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A mean-field model explaining most of the details in the magnetic phase diagram of ErNi2B2C is 
presented. The low-temperature magnetic properties are found to be dominated by the appearance 
of long-period commensurate structures. The stable structure at low temperatures and zero field is 
found to have a period of 40 layers along the a direction, and upon cooling it undergoes a first-order 
transition at ~ 2.3 K to a different 40-layered structure having a net ferromagnetic component 
of about 0.4 /ig/Er. The neutron-diffraction patterns predicted by the two 40-layered structures, 
above and below Tq, are in agreement with the observations of Choi et al. 
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A number of the (RE)Ni2B2C compounds are normal 
s-type superconductors with a of the order of 10 K. 
The superconductors are all of type II with a k in the 
range of 6-12. They are magnetic due to the rare-earth 
ions and in four of the compounds (RE = Dy, Ho, Er 
or Tm) the rare-earth ions are antiferromagnetically or- 
dered in the superconducting phase.EI The ordering wave 
vector Q in Er borocarbide is along an a axis and has 
a length of about 0.55 (in units of 27r/a), and the or- 
dered moments are along the a axis perpendicular to the 
ordering wave vector; the Neel temperature is = 6 
K and = 11 K. The rare-earth ions are placed in a 
body-centered tetragonal lattice, and in the case of Er 
a = b = 3.502 A and c = 10.558 A. It was proposed 
already in 1996 that the Er ions in ErNi2B2C develop a 
small ferromagnetic componervt in addition to the anti- 
ferromagnetic one below 2.3 kJ3 This makes the Er com- 
pound particularly interesting as presenting the case of 
a weak ferromagnetic state existing well below H^2- The 
ferromagnetic moment at 2 K was estimated to be 0.33 
jig per Er ion, which magnetization creates an internal 
magnetic field 47rM ~ 0.60 kOe close to the estimated 
value of the lower critical field H^^ . This opens up the 
possibility for the occurrence of exotic phenomena like a 
spontaneous vortex phase. Kawano et alu did detect a 
ferromagnetic moment below 2.3 K in a neutron diffrac- 
tion experiment, but saw nojsign of a spontaneous vortex 
phase. Recently, Choi et ala have made a detail neutron- 
diffraction experiment in which they measured all the 
higher harmonics of the antiferromagnetic structures oc- 
curring just above and below the Curie temperature, and 
they concluded that the structure at 1.3 K has a ferro- 
magnetic component of about 0.57 per Er ion. 

It is known that the superzone energy gaps on the 
Fermi surface induced by the antiferromagnetic ordering 
may have apStrong effect on the superconducting order 
parameter.l3i3 On the other hand, the superconducting 
order parameter does not seem to affect the antiferromag- 
netic ordering, but is of importance for a ferromagnetic 
system at low fields. The Anderson-Suhl mechanism,Ll 
that the Ruderman-Kittel-Kasuya-Yoshida (RKKY) of 
the rare-earth ions is strongly reduced in the long wave- 



length limit, has been demonstratedtl to be important 
in TmNi2B2C. Here, I shall mostly concentrate on un- 
derstanding the behavior of the magnetic moments in 
ErNi2B2C, and to start with I neglect the influence of the 
superconducting ordering on the magnetic properties. 

The crystal-field parameters of the Er ions have been 
determined by Gasser et al. from the crystal-field transi- 
tions observed by neutron scatteiing and from the high- 
temperature susceptibility data.l3 These parameters are 
used in the present work except that i?2 has been scaled 
by a small factor (1.08), see Table I. The ground state is 
a doublet, and an excited doublet is lying only about 0.6- 
0.7 meV above the ground state. This configuration leads 
to a four-clock behavior of the moments at low tempera- 
tures, i.e. the Er ions are easily magnetized along (100), 
they are hard to magnetize along the c direction, and 
when the field is applied along (110) the moments are 
(approximately) a factor y/2 smaller than the moments 
in the (100) case. 

Detlefs et al. have observed that the antiferromagnetic 
ordering in ErNi2B2C is accompanied by an orthorhom- 
bic distortion of the lattice, so that ajh — 1 — t^i ^ fu2 ~ 
2 • 10^^ at 3.7 K (when Q is along the a direction) .ETBe- 
cause of this observation, the following quadrupole cou- 
pling is included in the model: 

= -J2T. K [oTi^){OT) ~ UOT)'] (1) 

i m=±2 

The definition of the Stevens operators may be found in 
Ref. [|lOj. (O™) is the expectation value of the opera- 
tor averaged over all ions. The contribution to the free 
energy of the modulated quadrupolar moments (at the 
wave vector 2Q) is a factor of 100 smaller than that of 
the uniform term and is neglected. The value of is 



TABLE I: The Stevens operator parameters (meV). 

nO dO d4 dO d4 

i>2 i>4 i>4 i>6 

-0.0173 0.147-10"^ -3.3-10"^ -0.122 ■ lO"''^ 2.16 ■ lO"''^ 
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undetermined, but is only of minute importance for the 
model and is assumed to be equal K^. The calculated 
value of (O2) is 34.6 at 3.7 K. Assuming {c^ — Ci2)/2 
to be about 5 • 10^^ ergs/cm'^ (consistent with a De- 
bye temperature of the order of 400 K) then the dis- 
tortion observed by Detlefs et al. indicates a value of 
^ = {V/2N){c,, - Ci2)[(en - ^22)/ {Ol)f « 0.7 • lO'^ 
meV. This value is close to the one used in the final fit: 
= K-^ = 0.8 • 10-"^ meV. 
The two-ion interaction is assumed to be the sum of 
a Heisenberg interaction and the classical dipole-dipole 
interaction: 



3d Dal3iij)JtaJjl3 



(2) 



where the classical coupling is determined by 3^ — 
N{gfig)^ — 1.194 /ieV and the sum over the lattice of 



(3) 



The two-ion Hamiltonian is accounted for in the mean- 
field approximation, Jj-Jj- ^ Jj-(Jj)-l-(Jj)-Jj-(Jj)-(Jj). 
All the ordered structures are described by a wave vector 
Q along the a-axis and consist of ferromagnetic sheets 
perpendicular to Q. This means, firstly, that the different 
positions of the ions in the two sublattices have no direct 
consequences, corresponding to the use of a double-zone 
representation along (100), and, secondly, that only the 
total couplings between the different ferromagnetic layers 
are important in the model. The interplanar coupling 
parameters are 



E koj) + a 



D^l|.±(0j) 



(4) 



and the corresponding Fourier transforms 3|| ^(q)- The 
parameter 2 ± (q) denotes the coupling between the com- 
ponents of the moments which are lying in the a-h plane 
perpendicular to q. The case of involves the fol- 

lowing coupling parameters. 



= a±(LD) cos(0.558n7r) 



10,11,- •• ,16 



all defined in terms of 3^(LD). Notice, that the coupling 
parameter for n = 9 is included separately. 

The final values of the interplanar coupling parameters 
are shown in Table II, where d^{n) have been derived 
from the fitting of the experimental results discussed be- 
low, and 3|| (n)— have been calculated using Eq. (|^). 
The results for 3j^(q) and 3|| (q) are shown in Fig. 1. The 
most important parameters are the wave vector at which 
3j_(q) has its maximum, q = Qq = 0.558 a*, the maxi- 
mum value itseff 3j^(Qo) and aj^(O) = a||(0) = -12 ^eV. 
The value of the parameter 3_l(Qo) is assumed fixed at 
21.56 /J.eV, which leads to a mean-field value of the Neel 
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FIG. 1: The perpendicular and parallel components of the 
two- ion interaction in ErNi2B2C along [100]. 



temperature — 6.0 K. The classical interaction gives 
rise to a large difference tJ||(q) — 3j^(q) of the same or- 
der of magnitude as the Heisenberg term itself. Due to 
the singular behavior of the dipole sum this difference is 
cancelled at zero wave vector, i.e. 3|| (q) makes a jump by 
47r3^ = 15.0 /icV in the hmit of q ^ 0. 

The neutron-diffraction experiment of Choi et alu re- 
vealed a large third harmonic below 7^/2 showing that 
the modulation of the ordered moments approaches a 
square wave. The system shows a number of features 
indicating that commensurable effects are of decisive im- 
portance. The commensurable structures appearing may 
be rather complex limiting the possibility of deciphering 
the neutron-diffraction experiments. It has previously 
turned out to be of great value to assist the analysis 
of diffraction experiments by theoretical model calcula- 
tions of the stability of the structures which may occur. 
This has been the case in the study of the long-period 
commensurable atr.ur.tures in the elemental erbium and 
holmium metals.EElEiJ'Ej The method used here is the 
same as in the previous works just cited, i.e. the free ener- 
gies of different commensurable structures are calculated 
within the mean-field approximation by a straightfor- 
ward iteration procedure, and the results are compared 
to each other in order to identify the most stable struc- 
ture. 

The mean-field model presented above accounts for a 
great part of the observed properties of ErNi2B2C. Fig- 
ure 2 shows the magnetization curves calculated at 2 K in 
comparison with the experimental results. All the calcu- 
lated results are based on commensurable structures de- 
rived from the basic structure at Q = i (in units of a*). 
At low temperatures all the moments have a magnitude 
of about 7.9 /i^ and in one ferromagnetic layer they are 
either pointing parallel (u) or antiparallel (d) to the h axis 
(assuming Q || a). In the Q = \ structure the ferromag- 
netic layers perpendicular to the a axis are polarized sub- 
sequently uudduudd. . . along the a direction. Structures 
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TABLE II: The planar two-ion coupling parameters (/leV). 
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5.847 


-3.816 
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-0.650 
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-1.500 
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-14.286 


-3.024 


3.106 


-0.630 


0.250 


-0.084 


0.030 
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FIG. 2: The magnetization curves of ErNi2B2C at 2rM. The 
open circles are the experimental results of Cho et a/.tJ. The 
crosses connected by gashed lines show the experimental re- 
sults of Canfield et alB. The remaining solid and dashed lines 
are the calculated results. The experimental and theoretical 
results in the case of || [100] have been shifted upwards by 
2 units. 

with larger values of Q are derived from this structure 
by a periodic replacement of one or more of the uu (dd) 
double layers mith a single u (d) layer, the so-called spin- 
slip structures O Important structures in the present case 
are the Q — y[ structure consisting of the eleven lay- 
ered period d(uudd)^uu = d(5p), the eighteen layered 
Q — ^ structure with the period d(uudd)^u(dduu)^ = 
d(4p)u(4p), and the seven layered Q = | structure with 
the period duudduu = d(3p). Experimentally, the length 
of the ordering wave vector lies in the interval between 
0.554 and 0.548jj which range is covered by appropriate 
combinations of the three genetic structures. If the mo- 
ments are of constant length, the d(5p) and the d(3p) 
structures have a net ferromagnetic component equal to 
jj/^max and ^Mmax, respectively, whereas the structure 
d(4p)u(4p) has no uniform component. 

If the fundamental harmonic is the only one present 
the free energy has its minimum value at the wave vec- 
tor at which the two-ion coupling has its maximum, i.e. 
at Qo — 0.558 in the mean-field model defined above. 
However, the intensities of the higher-order odd harmon- 
ics increase rapidly as the temperature is lowered, which 
produce a shift of the ordering wave vector to a smaller 
value. In the model Q = 0.558 at T^, and it decreases 



rapidly between 5 and 3 K to become about 0.55 below 
3 K. This behavior is consistent with the experimental 
observationo of a change of the ordering wave vector be- 
tween 5 and 3 K from 0.554 to 0.548.. In the neutron- 
diffraction experiment of Choi et alu the peak due to 
the fundamental harmonic is found to be centered at 
Q = 0.548 at 2.4 K, however the higher harmonics indi- 
cate that the main part of the crystal is actually ordered 
in a commensurable structure with Q = = 0.55. At 
2.4 K, the stable structure with this wave vector is the 
40-layered structure d(4p)u(5p)u(4p)d(5p). This struc- 
ture, like the d(4p)u(4p) structure above, has no ferro- 
magnetic component. However, due to the increasing 
importance of the higher harmonics as the temperature 
is lowered these structures are found to become unstable. 
The present mean-field model predicts the occurrence of 
a first-order transition between the structures 

d(4p)u(5p)u(4p)d(5p) ^ d(3p)d(5p)d(5p)d(5p) (5) 

at 2.24 K. Numbering the layers from left to right by 1 to 
40, then the transition is accomplished by a reversal of 
the moments in the layers 9 and 20. The length of the mo- 
ments changes slightly from one layer to the next, and the 
ferromagnetic moment is calculated to be 0.33 /x^/Er just 
below the transition, and to be 0.40, 0.56 and 0.62 /i^/Er 
at 2, 1.3 and K^Jjespectively, in agreement with the ex- 
perimental valuegj'D of 0.33 /^g/Er at 2 K and 0.57 At^/Er 
at 1.3 K. The diffraction patterns of the different struc- 
tures have been calculated and are compared with the r& 
suits of the neutron-diffraction experiment of Choi et al.a 
in Fig. 3. At 1.3 K the d(3p)d(5p)d(5p)d(5p) configura- 
tion gives rise to both odd and even harmonics (the even 
ones are marked by arrows) twice as many as produced by 
the d(4p)u(5p)u(4p)d(5p) structure at 2.4 K. The only 
discrepancy of some importance is the large value calcu- 
lated for the intensity at /i = 0.7 at 1.3 K. The splitting 
of the peaks around h = 0.45 and h = 0.65 at 1.3 K may 
be explained if a minor part of the crystal is ordered in 
the Q = If structure, d(3p)d(5p)d(5p)d(5p)d(5p). No- 
tice, that even this small change of the fundamental Q 
(from 0.55 to 0.549) leads to easily observable modifica- 
tions in the positions of the higher harmonics. Hence, 
the overall agreement between the observed and calcu- 
lated positions of the higher harmonics indicate with a 
high degree of credibility that the main part of the crystal 
is a Q = structure at 1.3 and 2.4 K. There exist other 
choices for the low-temperature structure with Q = 
than the one proposed in (|^), but this structure, which 
is calculated to be the most stable one, is the only one 
producing a diffraction pattern that is reasonably similar 
to the one observed at 1.3 K. 
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FIG. 3: Scan along [h,0,0] at T = 1.3 K, 2.4 K, and 4.38 
K, measured with unpolarized neutrons. The data have been 
offset for clarity. The figure is a copy of Fig. 1 from Ref. in 
which the calculated diffraction patterns in the different cases 
(the heavy solid lines) have been incorporated. 

The model calculations indicate that 0(q) has a sharp 
peak at q = Qq, as constructed in terms of the parame- 
ter 0(LD). If this parameter is neglected the calculated 
variation of Q with temperature, or as a function of field, 
increases drastically. The dependence of Q on an applied 
field is goijig to be discussed in a forthcoming paper by 
Toft et alx3 In the present model all the structures be- 



tween Q = and | are so close in free energy below 3 
K that the model needs to be modified in order to dif- 
ferentiate clearly between the different Q values in this 
interval. This indicates that the peak in 3(q) is possibly 
even more pronounced than assumed in the present cal- 
culations. A strong peak in the RKKY-interaction may 
be produced by the nesting between differept areas on 
the Fermi surface discussed by Dugdale et aZ.J13 a nesting 
which is probably also responsible for the superconduct- 
ing properties of these compounds. 

It is important to realize that the ferromagnetic tran- 
sition in ErNi2B2C is not due to a ferromagnetic inter- 
action, which one is actually strongly negative. The 
ferromagnetic component is a byproduct of the low- 
temperature commensurable structure. Even a large 
change of 3(0) only have a slight influence on the transi- 
tion, as the exchange-energy gain of this phase relatively 
to the pure antiferromagnetic phase is determined by the 
combined contribution of all the even harmonics. This 
means that the influence of the superconducting electrons 
on this transitiorL_|as for instance through the Anderson- 
Suhl mechanismjSlj is weak. When the applied field in 
Fig. 2 is smaller than 1-2 kOe the response of the sys- 
tem is diamagnetic. The width of the fundamental har- 
monic in Fig. 3 indicates a correlation length along [1,0,0] 
which is at least 80 a, about twice the superconducting 
coherence length ^. The absence of a spontaneous vor- 
tex phase, i.e. the diamagnetic response of the system at 
low fields, therefore suggests that the magnetic correla- 
tion length perpendicular to Q in the Meissner phase is 
much shorter than the parallel component. The model 
is going to be examined more closely in a paper present- 
ing a neutron-diffraction detetmination of the magnetic 
structures in an applied field.113 

Valuable discussions with P. Hedegard, K. N0rgaard 
Toft, and N. H. Andersen are gratefully acknowledged. 
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